function transect_map(ctr,All);
transect=All.transect;
peak=All.peak;

%Load and plot grid data
R=region_data_all(ctr);
[lon, lat, z]=grdread2(['grdfiles/',R.name,'.grd']);
pl=find(lon>=R.cull_box(1) & lon<=R.cull_box(2));
lon=lon(pl); z=z(:,pl);
pl=find(lat>=R.cull_box(3) & lat<=R.cull_box(4));
lat=lat(pl); z=z(pl,:);

%load region data
fname=['../../generate_profiles/data/data_',R.name,'.mat']; 
load(fname,'C','P');

%plot map
h=surf(lon,lat,z,'FaceLighting','flat');
colormap(jet);
axis(R.cull_box);
if isempty(R.clim); R.clim=[min(z(:)) max(z(:))]; end;
caxis(R.clim);
shading flat;
for ct=1:3,
    light('Position',[-1 1 0]);
end;
h.AmbientStrength = 0.7;
h.DiffuseStrength = 0.7;
h=colorbar; h.Label.String='depth (meters)'; h.Label.FontSize=15; h.FontSize=15; 
font(gca,15); 
h=xlabel('longitude'); font(h,16);
h=ylabel('latitude'); font(h,16);
box on;

%plot transects
h4=axis;
ca=2;
for ct1=1:length(C.T) 
    %plot(C.T(ct1).ilon,C.T(ct1).ilat,'w','linewidth',4);
    text(C.T(ct1).ridge.rxy(1)+diff(h4(1:2))/50,C.T(ct1).ridge.rxy(2)-diff(h4(3:4))/50,num2str(ct1),'fontsize',20,'color','k');
    plot(C.T(ct1).ridge.rxy(1),C.T(ct1).ridge.rxy(2),'k.','markersize',36);
    
    %manually plot line spanning this study in black
    plbm=find(abs(C.T(ct1).t)<=780);
    plot(C.T(ct1).ilon(plbm),C.T(ct1).ilat(plbm),'k','linewidth',4);
    
    %plot locations of the B-M on each side where they exist
    plot(P.T(ct1).BM(1).lon,P.T(ct1).BM(1).lat,'m+','markersize',25,'LineWidth',4);
    if length(P.T(ct1).BM)>1
    plot(P.T(ct1).BM(2).lon,P.T(ct1).BM(2).lat,'m+','markersize',25,'LineWidth',4);
    end
end;
